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I.  INTRODUCTION 


High  speed  penetration  mechancs  involves  highly  complicated  physical 
phenomena,  all  operating  more  or  less  simultaneously,  depending  on  the  speed 
of  the  projectile  and  the  nature  of  the  target  material.  Classical  theories 
of  elasticity,  plasticity  and  visco-elasticity  are  often  Inadequate,  especially 
in  quantitative  terms.  In  explaining  these  penetration  problems.  Even  at 
moderate  speeds,  fracture, plastic  flow  and  new  surface  Initiation  occurs. 
Classical  discipines  are  not  designed  to  deal  with  these  types  of  physical 
phenomena . 

In  order  to  obtain  a  basic  understanding  of  penetration  phenomena, 
the  Isolation  of  fundamental  failure  mechanisms  and  their  rigorous  mathe¬ 
matical  treatment  are  necessary.  Moreover,  In  order  to  avoid  ad-hoc  assump¬ 
tions,  the  theoretical  underpinnings  of  such  an  approach  should  accommodate 
the  effects  of  microstructure  and  long  range  Interactions. 

To  achieve  the  above  goals,  the  nonlocal  continuum  theories  have 
been  developed  and  have  been  subjected  to  acid  tests  of  model  problems. 

The  results  have  been  most  gratifying.  After  much  laborious  efforts 
to  solve  complicated  integro-differentlal  equations  for  mixed  boundary- value 
problems.  It  has  been  possible  to  establish  simple,  meaningful  fracture  cri¬ 
teria.  Cricial  stresses  and/or  velocities  are  obtained  for  the  Initiation  of 
penetration  and  fracture.  The  estimates  for  the  cohesive  stresses  of  various 
materials  have  been  In  the  right  range,  I.e.  those  accepted  by  solid  state 
physicists  and  by  metal urgists. 

On  the  theoretical  side,  further  research  has  been  devoted  toward 
lR9>rov1ng  nonlocal  theories  and  fine  tuning  them  for  problem-solving  pur¬ 
poses.  The  correspondence  between  the  nonlocal  theories  and  solid  state 
physics  has  been  further  delineated.  Finally,  a  rigorous  thermodynamic 
derivation  of  nonlocal  plasticity  has  been  given. 
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Below,  first  we  present  certain  selected  research  results.  After 
brief  descriptions  of  the  solved  problems,  the  Important  conclusions  are 
restated.  In  the  next  section,  the  new  theoretical  work  In  nonlocal  formu¬ 
lation  Is  summarized. 

The  emerging  pattern  of  the  Interplay  of  the  theoretical  develop¬ 
ments  and  the  solutions  of  the  model  problems  Is  rather  encouraging,  in 
that  It  enables  us  to  formulate  a  unified  failure  criteria.  Further 
Intensified  research  Is  necessary  to  extend  these  basic  results  to  cases 
where  the  physics  of  the  penetration  phenomena  are  considered  In  more  detail. 


II.  RESEARCH  ACCOMPLISHED 

1.  NonLoaal  Fracture  Criterion 

A  plate  with  line  crack,  subject  to  uniform  tension  at  Infinity 
perpendicular  to  the  direction  of  line  crack  Is  a  crucial  problem  for  the 
understanding  of  the  rupture  of  brittle  materials.  The  classical  solution 
of  the  Griffith  problem  predicts  an  Infinite  hoop  stress  at  the  tip  of  the 
crack.  Because  of  this  singularity,  a  perfectly  adequate  criterion  of 
brittle  fracture  —  the  maximum  stress  hypothesis  —  used  for  all  other  types 
of  bodies  with  no  sharp  geometric  changes,  had  to  be  abandoned.  To  this  end 
Griffith  [1]  Introduced  his  celebrated  fracture  criterion  which  he  deduced 
by  equating  the  elastic  energy  to  "the  surface  tension  energy" 

f 

n)  -  C3  'g  ■ 

where  tQ  Is  the  applied  tension,  t  Is  the  half  crack  length.  The  Griffith 
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constant  Cq  Is  expressed  In  terms  of  Young's  modulus  E  ,  Poisson's  ratio 
V  *  and  the  surface  tension  energy  y  • 

The  Griffith  criterion  depends  on  the  surface  energy  concept  which 
itself  Is  not  Included  In  the  classical  elasticity  theory.  In  order  to  derive 
a  simple,  consistent  fracture  criterion,  a  better  description  of  the  crack  tip 
stresses  Is  necessary. 

Erlngen,  et  a1  [2]  solved  the  crack  tip  problem  by  means  of  nonlocal 
elasticity.  They  deduced  the  following  significant  results 

I)  The  stress  at  the  crack  tip  Is  finite  [Fig.  1].  It  Is  given  by 

(2)  tyy(£,0)  =  C(2£/a)^  tQ  C  =  0.73 

where  tQ  Is  the  applied  uniform  tension  at  Infinity. 

II)  When  the  maximum  stress  Is  equal  to  the  cohesive  stress  t^  ,  the 
fracture  will  occur.  In  this  case,  the  nonlocal  fracture  criterion  Is  given 
by 

(3)  t§£  *  (a/2c2)tj  -  Cg 

where  a  Is  the  atomic  distance.  Surprisingly  enough,  this  Is  the  Griffith 
criterion  with  the  added  benefit  that  the  Griffith  constant  Is  now  determined 
theoretically. 

III)  By  utilizing  the  formal  similarity  between  the  two  criteria  (1) 

and  (3),  we  can  predict  the  theoretical  cohesive  strengths  of  various  crystals. 
Once  the  experimental  value  for  y  Is  given,  the  calculations  can  be  carried 


out  easily  to  determine  the  cohesive  stresses.  In  Table  1,  the  nonlocal 
estimates  are  compared  to  those  of  atomistic  models  [3].  The  remarkably 
close  values  obtained  are  once  again  Indicative  of  the  far  reaching  power 
of  nonlocal  theories. 

2,  Initial  Penetraticm  by  a  Rigid  Pmoh 

One  of  the  most  fundamental  problems  of  penetration  mechanics  Is 
that  of  an  elastic  half  space  under  an  Indenting  punch  [Fig.  2].  The  clas¬ 
sical  solution  to  this  contact  problem  [4]  contains  stress  singularities  at 
the  edges  of  the  penetrator  no  matter  how  small  the  applied  loads  may  be. 
Hence  It  Is  very  difficult  to  predict  the  critical  stress  for  the  initiation  of 
penetration.  The  classical  elasticity  falls  to  apply  In  the  small  region 
around  the  edges  of  the  Indenting  block.  And  It  Is  exactly  this  region 
which  Is  crucial  to  our  understanding  of  the  penetration  phenomena. 

The  state  of  the  stress  around  the  edges  Is  heavily  Influenced  by 
the  microstructure  and  by  the  long  range  Interatomic  Interactions.  The 
recently  established  nonlocal  theories  Incorporate  these  effects  naturally 
[5].  To  this  end,  the  rigid  punch  problem  has  been  solved  by  means  of  the 
nonlocal  elasticity  theory.  The  stress  and  displacement  fields  are  obtained 
by  solving  a  set  of  complicated  dual  Integral  equations  [6]. 

The  variation  of  nonlocal  stress  at  the  surface  Is  plotted  In  Fig.  3. 
The  maximum  stress  Is  given  by 

I 

(4)  ^max^^O  “  ^  "  a/2kil 

where  tQ  Is  the  applied  average  stress,  k  Is  an  attenuation  factor  having 
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a  value  In  the  neighborhood  of  unity,  a  Is  the  atomic  distance  and  Zl 
Is  the  punch  width. 

The  nonlocal  analysis  of  the  rigid  punch  problem  leads  to  the 
following  conclusions: 

I)  Nonlocal  elasticity  predicts  finite  stresses  at  the  points  of 
contact  of  the  half  space  with  the  edges  of  the  penetrator. 

II)  By  eqcaiing  the  maximum  stress  to  the  cohesive  stress,  t^ 

(5)  p^^t  =  8  t^a/k  (k=l) 

It  has  been  possible  to  derive  a  penetration  Initiation  criterion  for  static 
loads  and  brittle  solids.  Note  that  In  arriving  at  this  result,  no  free 
parameter  adjustment  has  been  used. 

3.  Plug  Formation 

The  Impact  failure  of  a  target  plate  hit  by  a  projectile  occurs  by 
the  simultaneous  action  of  various  failure  mechanisms.  One  of  these  basic 
mechanisms  Is  that  of  plug  formation  (for  an  extensive  review  of  Impact 
phenomena,  see  e.g.  Backman  and  Goldsmith  [7]).  The  main  line  of  approach 
to  high  speed  penetration  problems  has  been  the  so-called  hydrodynamic  model 
which  employs  the  classical  viscous  fluid  theory  In  various  modified  and 
simpl'f^led  forms.  The  Inadequacy  of  classical  continuum  theories  In  des¬ 
cribing  the  high  gradient  phenomena  In  rupture  regions  has  already  been 
discussed. 


Thus  a  nonlocal  viscoelasticity  theory  has  been  Introduced  to  account 


for  the  high  strain  rate  effects  In  plug  formation  [8].  The  effect  of  the 
Impact  Is  represented  by  a  uniform  Initial  velocity  representation  over  a 
circular  region  of  the  plate  surface.  Since  the  plate  Is  very  thin,  only 
the  vertical  shearing  stresses  are  considered.  The  analytical  solution  to 
the  mixed  boundary  value  problem  yields  the  nonlocal  displacement  and  stress 
fields.  The  nonlocal  parameter  e  Is  determined  by  matching  the  Initial 
nonlocal  displacement  value  at  the  center  of  the  plug  (t  =  0,  r  =  0)  to  the 
experimental  result  given  by  Pytel  and  Davids  [9].  The  results  are  plotted 
In  Figures  4  and  5. 

The  time  decay  and  space  distribution  of  the  shear  stress  can  be 
used  in  design  considerations.  The  viscous  effects  cause  high  stresses  and 
small  displacements  In  the  first  few  microseconds  after  Impact.  As  the  time 
Increases,  the  maximum  stress  decreases  rapidly  and  the  displacement  continues 
to  Increase  at  a  slower  rate. 

The  present  solution  has  the  following  Important  characteristics: 

(I)  In  contrast  to  the  classical  analysis,  the  viscous  stresses  at  the 
periphery  of  the  plug  are  found  to  be  finite. 

(II)  The  critical  velocity  for  the  Initial  plug  formation  Is  given  by: 

(6)  »  (R/2.623  W)  tj. 

» 

where  t^  Is  the  cohesive  shear  strength,  u  Is  the  shear  modulus  and  R 
Is  the  projectile  radius.  Note  that  this  result  cannot  be  arrived  at  by 
the  classical  theory,  since  the  maxlmun  stress  Is  Infinite  and  therefore 
cannot  be  used  to  deduce  a  non-zero  Impact  velocity. 
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4.  Perforation  of  Viaoo-Plaatia  Platea 

The  high  speed  penetration  of  a  projectile  Into  a  plate  requires  the 
consideration  of  visco-plastic  effects.  In  addition,  the  nonlocal  Inter- 
molecular  forces  must  be  taken  Into  account.  In  [10],  a  nonlocal  analysis 
of  plug  formation  In  visco-plastic  plates  Is  given. 

The  model  problem  considers  an  Infinite  plate  of  thickness  hQ  being 
hit  by  a  cylindrical  rigid  projectile.  [Fig.  6].  The  Incidence  Is  normal 
and  the  Initial  velocity  Is  .  It  Is  supposed  that  after  such  an  Impact, 
a  cylindrical  portion  of  the  plate  with  radius  R  Is  set  into  uniform  mo¬ 
tion.  Since  the  plate  is  assumed  to  be  sufficiently  thin,  only  the  effects 
of  the  anti-plane  shear  stress  are  considered. 

The  nonlocal  analysis  generalized  the  classical  Bingham  solids.  Using 
appropriate  nonlocal  constitutive  equations,  which  include  the  plastic  yield¬ 
ing,  calculations  are  carried  out  to  obtain  the  stress  and  displacement  fields. 
The  Important  results  can  be  summarized  as 

(1)  The  maximum  value  of  the  shear  stress  Is  finite  and,  when 
It  reaches  the  cohesive  value,  the  penetration  begins. 

The  maximum  stress  is  given  by 

(7)  t^2  "  ■  ‘  |exp(-2/e)  I^(2/e) 

wnere 

(8)  e  =  a^/k^R^ 

and  I-j  Is  the  modified  Bessel's  function,  R  Is  the  radius  of 
the  projectile,  a  Is  the  atomic  distance,  k^^l  Is  the  attenua¬ 
tion  parameter  and  o  Is  the  ratio  of  the  yield  stress  to  the 
viscous  stress  (  uVq/R),  u  Is  the  shear  modulus. 
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(11)  Variations  of  plate  thickness  and  velocity  drop  with 

perforation  time  and  velocity  drops  with  plate  thickness 
are  plotted  In  Figs.  7  and  8.  Due  to  plastic  effects, 
longer  periods  are  required  for  complete  perforation.  The 
same  Is  valid  for  the  velocity  drop.  It  is  much  larger 
when  the  plastic  effects  are  considered. 


5.  Nonlocal  Elasticity  and  Lattice  Dynamics 

The  nonlocal  theory  of  elasticity  Incorporates  the  effects  of  the 
discreteness  of  the  lattice  and  of  the  long  range  atomic  interactions.  The 
similarity  between  the  underlying  physical  and  mathematical  structures  of 
the  lattice  dynamics  and  the  nonlocal  elasticity  can  be  very  useful.  For 
example,  the  continuum  formulation  of  nonlocal  theories  may  obliviate  the 
necessity  of  solving  rather  a  large  number  of  equations  of  motion  of  the 
discrete  lattice. 

The  correspondence  between  the  lattice  dynamic  and  the  nonlocal 
elasticity  becomes  especial ly  Important  for  materials  with  defect  structures. 
In  the  case  of  real  (imperfect)  materials,  the  equations  of  the  lattice 
dynamics  become  exceedingly  complex,  whereas  the  only  change  required  in 
nonlocal  field  equations  is  in  the  form  of  nonlocal  kernels. 

In  [11],  passage  is  made  from  nonlocal  elasticity  to  lattice  dynamics. 
Equations  of  motion  and  total  potential  energy  of  both  theories  are  compared 
and  found  to  be  Identical.  It  Is  further  demonstrated  that: 

(1)  The  nonlocal  elastic  moduli  can  be  obtained  by  matching  the 
dispersion  relations  of  plane  waves  In  nonlocal  elastic  solids 
to  those  of  lattice  dynamics. 


(11)  The  nonlocal  elasticity  theory  Intrinsically  Includes  the  surface 
effects  and  It  can  be  naturally  used  to  solve  problems  for  Imper* 
feet  finite  elastic  bodies. 


6.  Continuum  Meohanioa  at  the  Atomia  Scale 

In  [12]  and  [13],  the  difficulties  In  applying  the  classical  continuum 
theories  to  problems,  where  physical  or  geometrical  discontinuities  exist, 
are  delineated.  The  circumstances  when  the  nonlocal  theories  can  be  most  use 
ful  are  discussed. 

The  thermodynamic  restrictions  and  the  rigorous  development  of 
nonlinear  nonlocal  elasticity  Is  presented  with  special  emphasis  on  the 
functional  characteristics  of  the  constitutive  equations.  From  the  considera 
tions  given  there,  the  following  conclusions  can  be  derived: 

(1)  When  the  characteristic  length  of  the  external  disturbance 
\  Is  of  the  same  order  of  the  Internal  characteristic  length 
a  (e.g.  the  lattice  parameter),  the  classical  continuum  theories 
fail  to  apply.  One  of  their  basic  premises  X/a  >>  1  Is  no 
longer  true.  In  these  Instances,  If  one  Is  not  Interested  In  the 
Individual  motions  of  atoms,  nonlocal  theories  will  yield  suf¬ 
ficiently  accurate  results  for  the  cooperative  macroscopic  pheno¬ 
mena  of  the  real  bodies  (e.g.  cracks,  dislocations,  etc.). 

(11)  The  functional  character  of  nonlocal  constitutive  equations 

enables  one  to  distinguish  between  the  symmetry  of  the  molecular 
orientations  and  the  overall  symmetry  of  the  macroscopic  body. 

This  characteristic  of  nonlocal  theories  may  serve  as  a  powerful 
tool  In  cases  where  the  local  defect  syimetry  becomes  very 
Important  (e.g.,  bifurcation  of  running  cracks). 
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7.  Nonlocal  Plasticity 

The  nonlocal  character  of  plastic  deformations  and  of  plastic 
yielding  play  an  important  role  in  impact  and  penetration  problems.  The 
long  range  interactions  between  the  individual  dislocations  lead  to  the 
cooperative  phenomenon  of  yielding.  The  regional  loss  of  stability  and  the 
subsequent  plastic  flow  clearly  indicates  the  necessity  of  a  rigorous  thermo¬ 
dynamic  treatment  of  plasticity. 

In  [14],  a  theory  of  nonlocal  plasticity  is  developed  in  which  the 
state  of  stress  at  a  point  of  the  body  depends  on  elastic  ai^  plastic  strains 
at  all  points  of  the  body.  A  plastic  strain  measure  is  introduced  which  is 
not  necessarily  connected  with  the  external  strain  measure.  Thermodynamic 
treatment  given  in  [14]  generalizes  the  approach  of  Green  and  Naghdi  and 
their  coworkers  [15-17]  to  include  the  nonlocal  interactions.  Using  special 
types  of  functionals  (additive  functionals  in  the  sense  of  Friedman  and  Katz 
[18]),  constitutive  equations  are  obtained  for  anisotropic  and  isotropic 
materials. 

The  systematic  inclusion  of  nonlocal  effects  into  plasticity  theory 
provides  the  framework  for  the  mathematical  treatment  of  the  complicated  pro¬ 
blems  in  plasticity  and  microplastic  phenomena.  In  the  theoretical  formula¬ 
tion,  the  following  results  are  worth  noting: 

(i)  The  yield  criterion  is  derived  from  the  second  law  of 
thermodynamics.  The  "strain-space"  formulation  (i.e., 
using  the  elastic  and  plastic  strains  as  independent 
variables  rather  than  the  stress  and  plastic  strains), 
provides  a  self-consistent  set  of  constitutive 
response  functionals. 


n 


(11)  Constitutive  equations  are  obtained  for  the  rate  of  plastic 
strain. 

(Ill)  Special  constitutive  equations  are  given  for  the  nonlocal 
Isotropic  elastoplastic  solids,  with  no  work  hardening  and 
with  von  Mises-type  yield  criterion. 
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III.  TABLES  AND  FIGURES 


TABLE  I 

Cohesive  Stress.  4:  iJe~KY;KB8C’|t/«(I-v);«~interionicdisi.;  y  ■  Surface  energy. 


Eiperimental' 

Present  Theory 

Atomistic 

Models 

Type 

Crystal 

y 

CCS 

11*10-" 

CCS 

P 

A*' 

C 

K*10-" 

CCS 

Ky 

xlO-'* 

h 

KlO-" 

1,/E 

f,  £ 

fAl 

840 

2.51 

0.347 

2.86 

0.73 

5.216 

4.382 

1  2.38 

0  183 

rice  C- 

\Ni 

1725 

7.48 

0.276 

2.49 

073 

14.020 

24  185 

3.117 

0  163 

Body  C. 

Fe 

1975 

692 

0.291 

2.48 

0.73 

13.245 

26  158 

3.248 

0  182 

023 

Ionic 

LiF 

480 

4.40 

0.068 

2.014 

073 

6  407 

3.075 

1.2.36 

0  131 

Diam. 

C 

5400 

50.9 

0.187 

1.54 

0.73 

84.960 

458.782 

17.26t» 

0  143 

0.17 

Hei 

Zn 

575 

3.83 

0.333 

2.66 

0.73 

7.792 

4  481 

1.298 

0  1 2"’ 

Oil 

Atomistic  results  are  from  Table  7.1,  p.  160,  Lawn  and  Wilshaw  [3]. 


Figure  3:  Nondl Mens 1 one 1  Normal  Stress  at  the  Surface  [6] 
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Figure  8:  Velocity  Drop  Changes  with  Plate  Thickness  [10] 
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